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Abstract 

Oh' We calculate the wavelength dependence of the ratio of the linear polarization degree 

Q . to extinction (polarizing efficiency) P{A)/A(A) from the ultraviolet to near-infrared. The 

prolate and oblate particles with aspect ratios from a/b = 1.1 up to 10 are assumed to 
be rotating and partially aligned with the mechanism of paramagnetic relaxation (Davis- 
Greenstein). Size/shape/orientation effects are analyzed. It is found that the wavelength de- 
pendence of P(A)/A{A) is mainly determined by the particle composition and size whereas 
J> ' the values of P{A)/A{A) depend on the particle shape, degree and direction of alignment. 
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1 Introduction 



Modelling of the wavelength dependencies of interstellar extinction A(l) and linear 
polarization P{1) allows one to obtain information about such properties of inter- 
stellar grains as size, composition, shape, etc. Interstellar polarization also tells us 
about the structure of magnetic fields because it arises due to dichroic extinction 
of non-spherical grains aligned in large-scale Galactic magnetic fields [H], [El- A 
correlation between observed interstellar extinction and polarization shows that the 
same particles are responsible for both phenomena. 

Very often interstellar extinction and polarization are modelled separately (e.g., 0, 
BH). The modelling of these phenomena usually includes consideration of normal- 
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ized extinction E(A-Y) /E(B-Y) and normalized polarization given by Serkowski's 
curve (/'max is the maximum degree of polarization and /l^ax the wavelength corre- 
sponding to it, K is the coefficient) PiA)/P^^^ = exp[-K\n^{A^^^/ A)] (see, e.g., [51 
and discussion in [|6l, [|71). In these cases, the important observational data (abso- 
lute values of extinction and polarization) are ignored and normalized curves can 
be fitted using particles of diff"erent composition and slightly diff'erent sizes (see, 
e.g., discussion in yj, [|7]). For example, the consideration of absolute extinction 
gives the possibility to use the cosmic abundances to constrain grain models (see 
[[8|). In the case of interstellar polarization it is better to interpret the wavelength 
dependence of polarizing efficiency per unit extinction P(A)/A(A) (ratio of the lin- 
ear polarization degree to extinction). This allows one to estimate the degree and 
direction of grain alignment and particle shape. Note that previous consideration of 
polarizing efficiency was restricted by one wavelength band only (like P{Y)/A(Y) 
or P(K)M(K)). 

Here, we analyze the extinction and polarization produced by homogeneous spheroidal 
particles with diff'erent aspect ratios a/b (a and b are major and minor axes, respec- 
tively). Spheroidal particles have the important applications in various fields of sci- 
ence (spheroidal antennas, bacteria and microweeds, raining drops, etc.). By chang- 
ing a/b the particles with a shape varying from close to as sphere {a/b - I) up to 
needles (prolate spheroids) or disks (oblate spheroids) can be studied. The particles 
are assumed to be rotating and partially aligned with the Davis-Greenstein mech- 
anism (mechanism of paramagnetic relaxation) [9|. We focus on the behaviour of 
polarizing efficiencies and consider the normalized extinction A(/l) /Ay (i.e., ignore 
the cosmic abundances which will be taken into account in next paper). The the- 
ory is compared with observations of stars seen apparently through one interstellar 
cloud. 



2 Polarizing efficiency: observations 

Interstellar extinction grows with a radiation wavelength decrease while interstellar 
linear polarization has a maximum in the visual part of spectrum and declines at 
shorter and longer wavelengths (see [|71 and Fig. [H left panel). 

The polarizing efficiency of the diff"use interstellar medium is defined as the ratio 
of the percentage polarization (P) to the extinction (A) observed at the same wave- 
length P(A)/A(A). There exists an empirically found upper limit on this ratio 

/'maxM(V)^3%/mag, (1) 

where Pmax is the maximum degree of linear polarization which is reached on aver- 
age near the wavelength Aynax ~ 0.55 jum (see [|6l, [|71 for more discussion). 
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Fig. 1. Interstellar extinction (stars and thin line) and polarization curves for star HD 24263 
(left panel) and the polarizing efficiency of the interstellar medium in the direction of this 
star (right panel; solid curve shows the power-law approximation PjA oc A^-^"^). Observa- 
tional data were taken from [JO J (extinction) and 1.11,1 (polarization). 



We chose five stars located not far than 500 pc with measured ultraviolet (UV) 
polarization [11], then found the extinction (data published in [T0|| were mainly 
used) and calculated the ratio P{A)/A(A). It is suggested that these stars are seen 
through one interstellar cloud that is supported by a weak rotation of positional 
angle [fTTI . The obtained polarizing efficiencies are shown in Fig. |2l Apparently, 
first presentation of the observational data in the similar form was made by Whittet 
( [fT2l . Fig. 9) who gave the average normalized dependence P(A)/A(A) ■ Ay //'max- 
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Fig. 2. Polarizing efficiency of the interstellar medium in the direction of five stars. Obser- 
vational data were taken from ifTOl (extinction) and ifTTTl (UV and visual polarization), lfT3]| 
(IR polarization). 



Note that the polarizing efficiency generally increases with a wavelength growth if 
A I lum. It may be approximated using the power-law dependence P/A oc A"^. For 
stars presented in Fig. |2]the values of e vary from 1.41 for HD 197770 to 2.06 for 
HD 99264. 
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3 Modelling 



Let us consider non-polarized light passing through a cloud of rotating spheroidal 
grains. Interstellar grains are usually partially aligned (see, e.g., IfT?! ). Imperfect 
alignment is also supported by the fact that values of the polarizing efficiencies 
calculated for non-rotating or perfectly aligned particles are generally higher than 
the empirically estimated upper limit given by Eq. ([T]) IfTSJ . 

The extinction in stellar magnitudes and polarization in percentage can be found as 

m 

A(A) = 1.086A^d(Cext)i , PW = A^d<Cpoi)ilOO% , (2) 



where A^d is the dust grain column density and <Cext)^ and (Cpoi)^ are the extinction 
and linear polarization cross sections, respectively, averaged over the grain orienta- 
tions 

(C,,t)A = [l) jjj ^rlQ,,,f{^,/3)d^dcod/3, (3) 



7t/2 n njl 

<Cpoi)a = ^///^^'^ 2poi cos 2(A dipdojdp . (4) 





Here, ry is the radius of a sphere with the same volume as spheroidal grain, jS is 
the precession-cone angle for the angular momentum J which spins around the 
magnetic field direction B, ip the spin angle, a» the precession angle (see Fig. 1 
in m). The quantities Qext = (GSf + QlDl'^ and Q^,, = (2™ - Ql^^/l are, 
respectively, the extinction and polarization efficiency factors for the non-polarized 
incident radiation, /(^, ry) is the cone-angle distribution. 

We consider so-called imperfect Davis-Greenstein (IDG) orientation [9]. The IDG 
mechanism is described by the function f(^,/3) which depends on the alignment 
parameter ^ and the angle /3 

(^^cos^jS -I- sm y8)3/2 



The parameter ^ depends on the particle size ry, the imaginary part of the grain 
magnetic susceptibility (= HcOilT^, where is the angular velocity of grain), 
gas density n^, the strength of magnetic field B and dust {T^) and gas (Tg) temper- 
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2 ^ ^ where 6^"^ = 8.23 10^' n^yum. (6) 



The angle i/' in Eq. dH) is expressed via the angles (p, cjl),/3 and Q. (angle between the 
line of sight and the magnetic field, 0° < Q < 90°). Note that for the case of the 
perfect DG orientation (PDG, perfect rotational or 2D orientation) the major axis 
of a non- spherical particle always lies in the same plane. For PDG, integration in 
Eqs. ©, dH) is performed over the spin angle (p only. 

As usual, the particles of diff"erent sizes are considered. In this case the cross- 
sections (Cext)i and <Cpoi)i are obtained after the averaging over size distribution 
function n{rv) 

J (C)A(rv) n{ry) dry, 

where rymm and ry^ax are the lower and upper cutoffs, respectively. 
We choose the power-law size distribution function 

n{ry) oc Ty' , 

which was often used for the modelling of interstellar extinction [fT6ll . [fTTll . The 
average Galactic extinction curve can be reproduced using the mixture of carbona- 
ceous (graphite) and silicate particles (in almost equal proportions) with parame- 
ters: q = 3.5, rymin ~ 0.001 jum and rymax ~ 0.25 jum. 



4 Results and discussion 



Our calculations have been performed for prolate and oblate rotating spheroids 
of several sizes and aspect ratios. The particles consisting of astronomical silicate 
(astrosil) and amorphous carbon (ACl) were considered. The optical properties 
of spheroids were derived using the solution to the light scattering problem by 
the separation of variables method [fTHll and modem treatment of spheroidal wave 
functions [[T9l . 

Below we briefly discuss extinction by spheroids and in suflicient detail polarizing 
efficiency. The comparison with observations of two stars is performed in Sect. 14.31 
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4. 1 Extinction curve 



Extinction (and sometimes polarization) produced by spheroids was studied by 
Rogers and Martin (2011, Onaka [|2ll, Vaidya et al. B^, Voshchinnikov (231 and 
more recently by Gupta et al. [241. Mostly the particles with the aspect ratio a/b <2 
were considered. 

The parameters of our model (see Sectl3]) are: the type (prolate or oblate) and com- 
position of spheroids, the size ry or the parameters of the size distribution ry^mm, 
fv,max, 1, the degree specified through parameter Sq and direction of alignment Q. 
Figure [2] illustrates how variations of particle shape and alignment influence on the 
normalized extinction A{A)/Av. The value of for IDG orientation is typical of 
diffuse interstellar medium [fTSl . It is seen that the difference in curves are in evi- 
dence only in the UV. Prolate spheroids produce larger extinction when the angle 
of alignment reduces from ^l = 90° (magnetic field is perpendicular to the line of 
site) to Q = 0° (magnetic field is parallel to the line of site; Fig. [3l left panel). 
For oblate particles the dependence on Q. is opposite. However, only for PDG these 
changes are at a detectable rate. Extinction also grows when the particles (prolate 
and oblate) become more elongated or flattened (see Fig. [31 right panel). 




01234n678 
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Fig. 3. Normalized extinction for ensembles of aligned prolate spheroids from as- 
trosil. Spheroids have power-law size distribution with parameters: rymin - 0.001 //m, 
rv,max = 0.25 yum and q = 3.5; Sq = 0.196/im. The effect of variations of particle degree 
and direction of alignment (left panel) and particle shape (right panel) is illustrated. 

Variations of parameters of size distribution (rv^niin, fv,max, l) on extinction are sim- 
ilar to those for spheres are not discussed here (see, e.g., [25]). 



4.2 Polarizing efficiency: theory 

Theoretical dependence of the polarizing efficiency of the interstellar medium was 
discussed several times for particles of constant refractive index (see, e.g., Il20ll . 
(211 ). The wavelength dependence of P/A was considered in [|26ll for infinite and 
finite cylinders and in (Til for perfectly aligned spheroids of single size. Here, we 
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analyse firstly the dependence P(A)/A(A) for single size particles and then for en- 
sembles with a size distribution. 



4.2.1 Single size 

Some results for prolate particles are plotted in Figs. |4] and \5\ They show the po- 
larizing efficiency in the wavelength range from the infrared to far ultraviolet. The 
dependence P{A)/A{A) observed for two stars is given for comparison. We made 
calculations for compact grains and for porous grains using the Bruggeman mixing 
rule for refractive index and particles of same mass as compact ones. Note that cal- 
culated efficiencies can be considered as upper limits because some populations of 
grains (spherical, non-oriented) may contribute into extinction but not to polariza- 
tion. 





/ s 


astrosil, rv=0.1yLim 










prolate, a/b=2, 






IDG, «„=0.196Mm 








f 










0= 0° 

n=;io 






n=60° 






- ^ n=90° 











^1- -■„^', n=60° ' < 



\ 

\ 

- - - - s \ 


am. carbon, rv = 0.1,tim 
prolate, a/b-2, 
IDG, i5„=0.196Mm 


\ \ 




\ b 
\ \ 

O Q O 


n= 0° 

n=3o° 

n=8o 

-J ^ -^it ~ ~ ~ "=^''° 


o 




z 


4 6 8 1 


efficiency 


for homogeneous rotatinj 



spheroidal particles of astronomical silicate and amorphous carbon. The effect of varia- 
tions of particle composition and direction of alignment is illustrated. The open circles and 
squares show the observational data for stars HD 24263 and HD 99264, respectively. 

From Figs. H [5] one can conclude that the wavelength dependence of polarizing 
efficiency is mainly determined by the particle composition and size. Variations 
of other parameters influence on the value of efficiency (the dependence of P/A 
is scaled). A growth of the efficiencies P{A)/A{A) takes the place if we increase 
angle Q. (direction of alignment deviates widely from the line of site), parameter 
^0 (degree of alignment) or aspect ratio a/b (consider more elongated or flattened 
particles) and decrease particle porosity (volume fraction of vacuum) P or particle 
size ry. It is also evident that the contribution of particles with size ry « 0.1 //m 
to observed polarization must be rather important and the polarization produced by 
perfectly aligned particles significantly exceeds the observed one. 



4.2.2 Size distribution 

We also calculated the wavelength dependencies of P/A separately for spheroids 
from ACl and astrosil with standard power-law size distribution (rymin = 0.001 /urn, 
fv,max = 0-25 iJ.m and q = 3.5) and considered the effect of variations different model 
parameters. The results are shown in Figs.|6l and|71 Note again that the value of P/A 
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Fig. 5. Wavelength dependence of polarizing efficiency for homogeneous rotating 
spheroidal particles of astronomical silicate. The effect of variations of particle size, shape, 
porosity and degree of alignment is illustrated. The open circles and squares show the ob- 
servational data for stars HD 24263 and HD 99264, respectively. 

depends on the degree and direction of alignment and the shape of particles while 
the particle composition and size strongly influences on the wavelength dependence 
of polarizing efficiency. 
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Fig. 6. Polarizing efficiency ensembles of aligned prolate spheroids with ajh = 2 from as- 
trosil (left panel) and amorphous carbon (right panel). Spheroids have power-law size dis- 
tribution with parameters: rv,min = 0.001 jum, ry^max - 0.25 //m and q = 3.5; 6o = 0.196jum. 
The effect of variations of particle composition and direction of alignment is illustrated. The 
open circles and squares show the observational data for stars HD 24263 and HD 99264, 
respectively. 

As follows from these Figures, theoretical polarizing efficiency more or less resem- 
bles the observational dependencies only if the contribution of small particles into 
polarization is depressed (ry^nim > 0.05 /urn orq < 3, see Fig.|7]). It is well known for 
a long time that the polarization of forward transmitted radiation is determined by 
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Fig. 7. Polarizing efficiency ensembles of aligned prolate spheroids from astrosil. Spheroids 
have power-law size distribution, Q - 90° . The effect of variations of parameters of size 
distribution, degree of alignment and particle shape is illustrated. The open circles and 
squares show the observational data for stars HD 24263 and HD 99264, respectively. 

particles of small sizes while large particles produce no polarization independent 
of their shape (see discussion in [UJ and HZTlD . 

The influence of small particles on PjA is reduced if: a) the particles are absent in 
the interstellar cloud, b) the small particles are less oriented or c) the shape of small 
particles is closer to spherical. The reason a) must manifests itself by the depressed 
UV extinction while the reasons b) and c) would be difficult to distinguish. 

The relationships found can be used for more detailed comparison of the theory 
with observations. 



4.3 Comparison with observations 



We calculated the normalized extinction A(/l) /Ay and the polarizing efficiency and 
compared the results with observations of stars mentioned in Sect. [21 In order to 
fit the UV bump on the extinction curve a small amount of graphite spheres with 
radius ry ^ 0.015 jjm was added to the mixture of silicate and carbon spheroids. 

According to the model presented in Sect. [31 we adopted that the same particles 
produce extinction and polarization and the alignment of silicate and carbon parti- 
cles is identical. At the same time, parameters of size distribution for silicate and 
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carbon can differ. 



The results of comparison are shown for two stars in Figs.[8]and|9l 
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Fig. 8. Normalized extinction (left panel) and the polarizing efficiency (right panel) in 
the direction of the star HD 147933. The error bars show Icr uncertainty. Dashed curves 
show the results of calculations for model with prolate spheroids with a/b = 5 from 
amorphous carbon (24%) and astrosil (73%) and graphite spheres (3%). Spheroids have 
a power-law size distribution with the parameters: amorphous carbon, rymin = 0.03 jum, 
fv,max — 0.15//m, q = 1.5; astrosil, fv^min — 

0.08 /^m, fy^max ~ 0.25 jum, q — 1.5. Alignment 

parameters are: 60 = 0.5yum, Q. = 35°. 

The coincidence of the theory with observations is rather good. The better agree- 
ment can be achieved if we assume that the alignment parameters for silicate and 
carbonaceous grains differ. For example, it is possible to consider non-aligned car- 
bonaceous grains (as it was made in ^2M ). Evidently, more attention should be 
given to modem alignment mechanisms. 
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Fig. 9. Normalized extinction (left panel) and the polarizing efficiency (right panel) in the 
direction of the star HD 62542. The error bars show Icr uncertainty. Dashed curves show the 
results of calculations for model with prolate spheroids with afh = 2 from amorphous car- 
bon (85%) and astrosil (10%) and graphite spheres (5%). Spheroids have a power-law size 
distribution with the parameters: amorphous carbon, rymin - 0.07 /im, rymax = 0.10//m, 
q = 2.7; astrosil, rymin - 0.005 //m, rymax = 0.25 //m, q = 2.7. Alignment parameters are: 
5o = 0.5//m, a-40°. 

These improvements of the model require to invoke additional observational data 
such as dust-phase abundances in the directions of considered stars, information 
about magnetic fields, etc. Otherwise, the modelling cannot give non-ambiguous 
results. In particular, the important feature is the wavelength where the ratio P/A 
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reaches a maximum (see Figs. In order to find this quantity we need the 

observations of polarization in the near-IR (J, H, K bands) which are absent for 
some considered stars (Fig.O. 



5 Conclusions 

In the frame of the model of partially aligned spheroidal grains we studied the 
wavelength dependence of the ratio of the linear polarization degree to extinction 
(polarizing efficiency) P(A)/A(A). 

The main results can be summarized as follows. 

1. It is found that the wavelength dependence of P{A)/A{A) is mainly determined 
by the particle composition and size. 

2. The values of P{A)/A{A) depend on the particle shape, degree and direction of 
alignment. 

3. The modelling of the only wavelength dependence of polarizing efficiency does 
not allow one to determine all parameters dust ensemble. Therefore, interpretation 
of the observations must include consideration of extinction. 
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